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Abstract 

The Hamiltonian formulation plays the essential role in constructing the framework of modern 
physics. In this paper, a new form of canonical equations of Hamilton with the complete symmetry 
is obtained, which are valid not only for the first-order differential system, but also for the second- 
order differential system. The conventional form of the canonical equations without the symmetry 
[Goldstein et al, Classical Mechanics, 3rd ed, Addison- Wesley, 2001] are only for the second-order 
differential system. It is pointed out for the first time that the number of the canonical equations 
for the first-order differential system is half of that for the second-order differential system. The 



Or 



nonlinear Schrodinger equation, a universal first-order differential system, can be expressed with 
the new canonical equations in a consistent way. 
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The Hamiltonian viewpoint provides a framework for theoretical extensions in many 
areas of physics. In classical mechanics it forms the basis for further developments, such as 
Hamilton- Jacobi theory, perturbation approaches and chaos [1]. The canonical equations of 
Hamilton in classical mechanics are expressed as [l| 

1i = g^7> (i= l,-- - ,n), (1) 
-Pi = q^, (« = !,•••,"), (2) 

where g« and Pi are said to be the generalized coordinate and the generalized momentum, 
qi = dqi/dt is the generalized velocity, pi = dpi/dt, and H is the Hamiltonian of the system. 
The generalized momentum pi is defined as pi = Jjp where L is the Lagrangian of the system. 
And the Hamiltonian H is obtained by the Legendre transformation H = Yl7=i QiP* ~ L. A 
set of (qi,Pi) forms a 2n-dimensional phase space. 

The canonical equations ([1]) and (|2J) are not only valid for the discrete system, but also 
can be extended to the continuous system as [j| 

Qs = ^, (s = l,---,N), (3) 

OPs 

-Ps = ^, (s = l,---,N), (4) 
oq s 

where the subscript s represents the components of the quantity of the continuous system 
hereafter it all denotes 1, • • • ,N, |^ = |^ - |-^- and f- = |^ - denote the 



functional derivatives of h with respect to q s and p s with q SjX = |p and p StX = j^-, q s and 



Sq s dq s dx dq s ,x Sp s dp s dx dp s 

~ a and j)q t ■ — 

OX ca,j, g x 

p s are the generalized coordinate and the generalized momentum, respectively, and h is the 
Hamiltonian density of the continuous system. Like the discrete system, the generalized 
momentum p s for the continuous system is defined as 

Ps = TTT - ) (5) 

dq s 

and the Hamiltonian density h for the continuous system is obtained by the Legendre trans- 
formation as 

N 

h = ^2 QsPs - I, (6) 



where I is the Lagrangian density. But it is significantly different for the continuous system 
that q s and p s are now not only functions of time t, but also the spatial coordinate x, where 



2 



the spatial coordinate x is not the generalized coordinate, but only serves as the continuous 
index replacing the discrete i in Eqs.flTJ and (j2J). To distinguish t from the coordinate x, 
we refer time t as the evolution coordinate. p s and q s define the infinite-dimensional phase 
space, h is a function of q s ,p s and q SjX but not p s>x , so = J^-, then Eq.flH]) can be also 
expressed as 

dh 

= J5~- 7 
"The advantages of the Hamiltonian formulation lie not in its use as a calculational 
tool," as has been pointed out in ref.jl], "but rather in the deeper insight it affords into the 
formal structure of mechanics." We have a greater freedom to select the physical quantities 
to designate as "coordinates" and "momenta", which have equal status as independent 
variables. As a result we can present the physical content of mechanics in a newer and more 
abstract way. The more abstract formulation is primarily of interest to us today because 
of its essential role in constructing the framework of modern physics, such as statistical 
mechanics and quantum theory. 



Canonical equations of Hamilton with symmetry 

The Newton's second law of motion in classical mechanics, based on which the Hamilto- 
nian formulation is established, is the second-order differential equation about the evolution 
coordinate (here the evolution coordinate is time), which is in this paper referred as the 
second-order differential system. Similarly, the first-order differential system is the first- 
order differential equation about the evolution coordinate. The Lagrangian density of the 
second-order differential system of the continuous systems is expressed in general as [l| 

N N N 
s=l k=l 8=1 

where A s k, B s , C depend on not only q s but also q SjX in general. The generalized momentum 
can be obtained by the definition (j^D as 

N 

p s = ^(A sfc + A ks )q k + B s , (9) 
fe=i 

which is a function of q s , q s and q S)X . The number of Eqs. © is N, and there are 4iV variables, 
which are q s , q s , p s and q StX . So the degree of freedom of Eqs. (jUJ) is 3iV. Then q s ,p s and 



q S:X are taken as independent variables. The generalized velocities q s can be expressed with 
these independent variables. 

Besides the second-order differential systems, there are a number of the first-order dif- 
ferential systems to model physical phenomena, for example, the nonlinear Schrodinger 
equation (NLSE) 

dip ld 2 ip , 9 , , 

which is a universal model that describes many nonlinear physical systems and can be 



5], nonlinear acoustics |6j, Bose-Einstein 
5|, the NLSE ( ITU]) governs the propagation 



applied to hydrodynamics |2|, nonlinear optics 
condensates and so on. In nonlinear optics {3 
of the slowly-varying light-envelope, where the evolution coordinate t is the propagation 
direction coordinate, the light-envelope if is a cw paraxial beam in a planar waveguide {^J or 
a narrow spectral-width pulse in optical fibers {4], 0], and x is a transverse space coordinate 
for the beam and a frame moving at the group velocity (the so-called retarded frame) for 
the pulse, respectively. We will show in the paper that the canonical equations for the 
first-order differential system are significantly different from Eqs.(j3]) and which are the 
canonical equations only for the second-order differential system. But in ref. l|, the canonical 
equations (EJ) and (@| were obtained regardless whether the continuous system is the first- 
order differential system or the second-order differential system. In ref. (4], the canonical 
equations for the NLSE were considered to be the same as Eqs.([3]) and (jJJ). But obviously 
the canonical equations (|3]) and (j4]) are not valid for the NLSE, from which it is impossible 
to obtain the NLSE or its complex-conjugate equation, as will be shown in the following. 
Attempt was made to deal with the difficulty in ref.js], but the canonical equations for the 
NLSE they obtained were inconsistent, as will be shown in the following. In this paper we 
obtain the canonical equations with the complete symmetry, which are valid not only for 
the first-order differential system, but also for the second-order differential system. 

For the first-order differential system of the continuous systems, the Lagrangian density 
must be the linear function of the generalized velocities q s . If the Lagrangian density is a 
quadratic function of the generalized velocities like Eq.flB]), the equation of motion, i.e., the 
Euler-Lagrange equation 

»«--«. = (11) 
dt dq s 5q s 

will be the second-order differential equation about the evolution coordinate t, which is 
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in contradiction with the definition of the first-order differential system. Therefore, the 
Lagrangian density of the first-order differential system can only be expressed as 



Besides, R s in Eq. (jl2j) is not the function of a set of q SjX . If R s is also the function of a set 
of q S;X , there will be such terms as q s>x q s appearing in Eq. f[T2"j) . Substitution of Eq. f[T2"j) into 
Eq. fllll) leads to the appearance of the mixed partial derivative term the term can be 
changed to -^f by the rotation of the coordinate frame. After the coordinate transformation, 
the Euler-Lagrange equation fflTl) is changed to the second-order partial differential equation 
of the standard form. Then the system expressed with the standard form is a second-order 
differential system. Consequently, the generalized momentum p s , which is obtained by the 
definition (j^D as 



is only a function of q s . This is of significant difference from the case of the second-order 
differential system, where the generalized momentum p s is the function of not only q s , but 
also q s and q StX , as shown in Eq.(|9]). There are 2N variables, q s and p s , in Eqs. (fTBl . And 
the number of Eqs. ( {TBI is N, which also means there exist N constraints between q s and 
p s . So the degree of freedom of the system given by Eqs. (TT31 is N. Without loss of 
generality, we take q±, • • • ,q u and Pi, • • • , p M as the independent variables, where u + n = N. 
The remaining generalized coordinates and generalized momenta can be expressed with 
these independent variables as q a = f a (qi, ■ ■ ■ , q v ,Pi, ■ • • ,P^)( a = u + 1, - • - , N), and pp = 

We now derive the canonical equations for the first-order differential system. The total 
differential of the Hamiltonian density h can be obtained by using Eq.fjH]) 



N 




(12) 



8=1 



p s = R s (q s ), 



(13) 




N 



(14) 



while the total differential of the Lagrangian density l{q s ,q s ,q s ,x) with respect to its argu- 
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ments is 



A=l 



^9/ * dl 

+ 5^ s+ 5^^- (15) 



Substitution of Eq.(|I5J) into Eq.dUD yields 



( ■ d 9P dl df a dl\ 
dh=y > gp- > — — \dq x 

+EK + E ^ff- E &w) dPv 

N ^ 

"S^" 9 " (16) 

Since the total differential of h(qi, ■ ■ ■ , q v ,Pi, 4 4 4 ,P/i, <?s,z) with respect to its arguments can 
be written as dh = J2\=i |^<7a + EJ=i J^^ + ZlfLi -^dq SjX , by comparing this equation 
with Eq. (ll6jl . we obtain 2N equations 



dh ^ . dgp A dl df a dl 

^r^wr^wrw^ (17) 
^+E^-E^ as) 

dh dl 

(19) 



dq s , x dq S:X ' 

where A = 1, • • • , v, T] = 1, • • • , /i, s = 1, • • • , JV. From the Euler-Lagrange equations ( |TTi) . 
we obtain 

dq s dt dq s dx dq s , x ' 
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Substituting Eq.flU]) into Eqs.JIZD and (HD, we have 

dq\ ^ X ^ dq\ ^ ^° dq\ 

y 19 d/ gj« g 91 (21) 

dfr _ . ■ ^ _ y ■ d I^. 

y- d__dl_dfa , 22 , 
a t? +1 dx dq a , x d Pri ' 

Then substituting Eq. (TT^|) into Eqs. (l2Tj) and (1221) . we obtain the N canonical equations for 
the first-order differential system 

Sh_ = _. ^ . c*9i _ . dfa 

+ y (23 ) 

^-f dx dq a , x dq x ' 
Sh _. ^ . dgj3_ _ . dfa 

+ y d_^d^_ (24) 

^ dx dq a , x dpr, 

To obtain Eq. (1241) . we have used 4^- = J^, because h is not a function of p v>x . The canonical 



equations ( 1231) and ( |24|) can be expressed in a completely symmetric form as 

Sh sr^ f ■ 9 Ps ■ 9q s \ v-^ d dh df a 

ras " Ps ^ j£, SaS^' (25) 

—=y ( ■ 9ps 9q A + d__dh_df^ , 2Q , 

(A = I, - •• ,v, rj = I, - •• and i/ + /x = N), because p x = Ea'=i PA'I^S = S?'=i^'^"' 
X^=i = 0' and Xa=iPA§^ = 0- The canonical equations ( 125]) and (|26|) can be easily 
extended to the discrete system, which can be expressed as 

dh -A / . dp s . dq s \ 

wrh{ q 'wr p 'wJ- (27) 

dh ^ f . dp s . % \ 

a^M^-^J' (28) 
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where A = 1, • • • , u, rj = 1, • • • , fi, and u + n = N . 

Although they are derived from the first-order differential system, the canonical equations 
of Hamilton (125]) and with symmetry can be proved to be reduced to Eqs.([H]) and (0J 
without symmetry, if all the N generalized coordinates q s and N generalized momenta p s 
in Eqs.([25D and <$Mj are independent, because |^ = |^ = ^ = gs. = 0(a = z/ + 
1, • • • ,N, (3 = fi + 1, • • • , N) in this case. In fact, this is just the case of the second- 
order differential system. So the canonical equations fl25|) and f )26|) with the completely 
symmetry are valid not only for the first-order differential system, but also for the second- 
order differential system. The conventional form of the canonical equations ([3]) and (jlj) 
without the symmetry are only for the second-order differential system. In addition, the 
number of Eqs. f[25"j) and (126]) is N in the case for the first-order differential system, and is 
half of that of Eqs.(j3]) and fll]). This can be explained in the following way. For the second- 
order differential system, the Euler-Lagrange equations ( ITT]) are the second-order differential 
equations about the evolution coordinate, where the Lagrangian density is expressed as 
Eq.flHJ). The number of Eqs.f llip is N. It is well known that one second-order differential 
equation can be reduced to two first-order differential equations {9]. Then 2N canonical 
equations, which are the first-order differential equations about the evolution coordinate, 
can be obtained from the N Euler-Lagrange equations. But it is significantly different for 
the first-order differential system that the Euler-Lagrange equations (|TT|) are the first-order 
differential equations about the evolution coordinate, so only iV canonical equations are 
obtained from Eqs.(|TTi). 



Hamiltonian formulation for the NLSE 



After constructing the canonical equations of Hamilton for the first-order differential 
system, we discuss its application to the NLSE. It is known that the Lagrangian density for 
the NLSE can be expressed as Q I = -|(^*^ - + §l§: | 2 - |M 4 - The NLSE is 

complex, and therefore it is an equation with two real functions, the real part of ip and its 
imaginary part. It is convenient to consider instead the fields ip and <p>* which are treated 
as independent from each other. Therefore, iV (the components of the quantity) for the 
NLSE equals two, i.e., there are two generalized coordinates, q\ = ip* and q 2 = ip, and two 
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generalized momenta can be obtained by using the definition (j5J as 



Pi = = ( 29 ) 



The Hamiltonian density for the NLSE can be obtained by use of Eq.flB]) as [4] 



h = -i 

2 



dip 



dx 



+ (30) 



If the generalized coordinate q\ and the generalized momentum pi are taken as the inde- 
pendent variables, the remaining generalized coordinate q 2 and the remaining generalized 
momentum p 2 can be expressed via the relations fl29|) as q 2 = —2ipi and p 2 = — respec- 
tively. We should also note that the Hamiltonian density h is also the function of q SiX , which 
are independent from q s and p s . Then for the NLSE, the Hamiltonian density ( 130]) should 
be expressed with the independent variables qi,Pi,qi, x and q 2tX as 

h = -~qi, x q2,x - ZQiPi- (31) 

We should note that v — p, — 1 and N = 2 for the NLSE, which means that the equations 
(125]) have only one equation, so do Eqs.()26l). Therefore, the canonical equations (1251) and 
(J2SJ) will yield two equations for the NLSE. The left side of Eq.( 125|) is obtained as 



5h _ 1 d 2 cp 
5qi 2 dx 2 

and its right side is 



+ M V (32) 



. , . dp 2 . . 
-Pi + ?2^— = (33) 

Then the NLSE ( FlOl) can be obtained. Using the other canonical equation, the left side of 
Eq. fl26l) is obtained as 

^- = -Mlpi = -2iM V, (34) 

and its right side is 

. dq 2 d dh dq 2 .d 2 ip* 
qi - p Wi + dx-dq^dp-i =2V (35) 
Then the complex conjugate of the NLSE is also obtained. Therefore, the canonical equations 

( 125]) and ( 12 6p are consistent in the sense that the NLSE can be expressed with one of the 

two canonical equations, and its complex conjugate can be expressed with the other. 

In ref.[4], the canonical equations for the NLSE were considered to be the same as those 

for the second-order differential system. Then, according to the definition (j5J), p v = 31/ 'dip, 



p v must be —i/2(p* but not — iip*. It was artificially doubled in ref. {4] so that p v = —i(p* 
in Eq. [5.1.29] (to avoid confusion, we replace the parentheses by the brackets to represent 
the formulas in the references) to make the NLSE derived from the canonical equation (J3J). 
In fact, substitution of the Hamiltonian density f[3"Ul) into Eq.(HJ) only yields + |V5_y? + 
\<p\ 2 <p = 0, which in fact does not be the NLSE ( ITUl) . In ref. 8|, the canonical equations the 
authors obtained are Eqs.[3.87] and [3.81], the latter is the same as Eq.(J7]). Although the 
NLSE can be derived from Eq.[3.87], its complex conjugation could not be obtained from the 
other, Eq. [3.81]. We now show this claim. Substituting the Hamiltonian density (130!) into 
Eq.©, the left side of it is obtained as ^jf , and the right side is ^ = = -2i\<p\ V , 



where p v * = Then the equation — | -ff +M = can be obtained, which is absolutely 
not the complex conjugate of the NLSE. Therefore, the canonical equations for the NLSE 



obtained in ref. 



Conclusion 
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are inconsistent. 



Summarizing, we obtain a new form of canonical equations of Hamilton with the complete 
symmetry, which are valid not only for the first-order differential system, but also for the 
second-order differential system. The conventional form of the canonical equations without 
the symmetry are only for the second-order differential system. The number of the canonical 
equations for the first-order differential system is N rather than 2N like the case for the 
second-order differential system. The canonical equations for the NLSE are two equations, 
which are consistent in the sense that the NLSE can be expressed with one of the canonical 
equations and its complex conjugate can be expressed with the other. The Hamiltonian 
formulation for the continuous system can also be extended to the discrete system. 



It is well known that the symmetry plays an important role in theoretical physics 11 ]. 
The search for and the discovery of new symmetries promote the exploration of fundamental 
laws of physics. Based on the idea, the conical equations of Hamilton with the complete 
symmetry found by us might find their appropriate position in modern theoretical physics. 
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